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Q . We ensure successful explosions (of otherwise non-explosive models) by enhancing the 
neutrino luminosity via reducing the neutrino scattering cross sections or by increasing 
^ I the heating efficiency via enhancing the neutrino absorption cross sections in the heating 
region. Our investigations show that the resulting electron fraction Yg in the innermost 
ejecta is close to 0.5, in some areas even exceeding 0.5. We present the effects of the 
resulting values for Ye on the nucleosynthesis yields of the innermost zones of core collapse 
supernovae. 



1. Introduction 

Presently, supernova explosions cannot be explained in ID spherically [ 121 EI or 2D 
rotationally [ H) |3] symmetric radiation-hydrodynamics calculations. Besides the need 
for further refinement and consolidation of the current simulations, successful explosion 
simulations may rely on rotation and magnetic fields [ E] or unknown microphysics, 
e.g. for the neutrino interactions. For the prediction of the nucleosynthesis yields in 
supernovae, we face the problem how to simulate the behavior of ejected matter real- 
istically, given the existing problems to obtain explosions in self-consistent models with 
accurate neutrino physics. Artificially induced explosions for supernova nucleosynthesis 
predictions [ |H1 El HD] are a valid approach for the outer layers, but are inconsistent for 
the innermost layers, affecting the Fe-group composition. The nucleosynthesis in these 
zones (undergoing Si-burning) is dominated by the electron fraction Yf.. The constraints 
from observations require the resulting Y^ in the innermost zones to be > 0.5 [E]. 
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Figure 1. (a) The upper panel shows the time evolution of the electron fraction in a fluid 
element close to the mass cut (thick solid line). The dotted line is the equilibrium electron 
fraction for infinitely long exposure to the neutrino field. The dashed and thin solid lines 
represent equilibrium with only neutrino absorption or emission considered. The electron 
fraction is determined by competition between the neutrino interaction rates and the 
matter ejection timescale. (b) Electron fraction as function of mass when the runs are 
stopped (0.4s - 0.6s after bounce) with (A) reduction of the neutral current scattering 
opacities by a given percentage (solid lines), or (B) enhancement of the reaction timescales 
(dashed lines). The open circles are final electron fractions for model A(40%). 



2. Hydro dynamical Simulations 

The simulations are performed in spherical symmetry with general relativistic Boltz- 
mann neutrino transport using the code AGILE-BOLTZTRAN. For a detailed description 
of the code refer to [IHj. We use two different approaches to simulate explosions. In 
the first approach, the neutral current scattering opacities of neutrinos on free nucleons 
are scaled with a factor ranging from 0.1 to 0.7 while using standard cross sections for 
electron/positron/neutrino/anti-neutrino capture on free nucleons (series A ^). This leads 
to a faster deleptonization of the protoneutron star such that the neutrino luminosities 
are boosted in the heating region, as it might happen with rigorous protoneutron star 
convection or dramatic changes in the input physics at high densities. In the second ap- 
proach, the charged current neutrino reactions (and their inverse reactions) are increased 
by equal factors in the heating region (series B) to crudly mimic a prolonged exposure to 
neutrino heating in the convective heating region which cannot be appropriately treated 
in spherically symmetric simulations. This enhances the time scale for neutrino heating 
without changing the equilibrium Ye. 

The models are based on a progenitor of 2OM0 [E]. For each series the choice of 



^Note that series A has been calculated at reduced neutrino angle resolution which additionally enhanced 
the neutrino abundances. 
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Figure 2. Mass fractions of the stable nuclei for model A(40%) normalized to solar abun- 
dances. For both calculations shown the postprocessing is based on the same hydro profile. 
The differences originate from the set of reactions included in the nuclear reaction net- 
work: no neutrino-induced reactions included, or neutrino captures on free nucleons and 
on nuclei included. 



parameters is such that they correspond to a barely exploding model, an extremely ex- 
ploding model, and a model with average explosion properties. The explosion energies 
are of the order of lO^^ergs. After ~ 100ms the accretion front stalls in all models and 
the explosion is launched after a delay of additional ~ 100ms. Three different phases can 
be distinguished in the mass trajectory: During the first phase the mass element falls 
inwards in the gravitational potential. This phase lasts ~ 200ms and ends with an almost 
instantaneous deceleration as it falls through the accretion shock. For approximately 
200ms it drifts around in the heating region until, in a third phase, it is ejected to larger 
radii lifting the electron degeneracy. When the electron chemical potential falls below the 
mass difference between neutrons and protons, protons are favored over neutrons [ . 
Figure ^ compares the evolution of the electron fraction to equilibrium values obtained 
under different assumptions. The simulations are stopped when convergence problems 
appear due to the mass separation between the remnant and the ejecta. 

3. Nucleosynthesis 

For the nucleosynthesis results presented here, we consider only the first few zones out- 
side of the mass cut. The position of the mass cut emerges consistently from the simulation 
as the region of bifurcation in which the density has dropped below ~ lO^g/ cm^. Based on 
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the hydrodynamics profiles, the detailed nucleosynthesis is calculated in a framework of 
consistently treating all weak interactions, namely neutrino/anti- neutrino absorption on 
free nucleons, the inverse reactions thereof [electron/positron capture], and neutrino/anti- 
neutrino absorption on nuclei. 

All of the models develop a proton-rich region around the mass cut. This can be seen in 
Figure^b) which shows the electron fraction ~ 0.5s after bounce. The open circles denote 
the final electron fraction of the exemplary run A(40%) at temperatures T < 0.1 x lO^K. 

Figure 121 shows the mass fractions of the stable nuclei for model A(40%). In this model, 
the mass cut sits at l.SllM©. For the other models that mass cut is between 1.44M0 and 
1.586M0. The light elements are mainly produced in the outer zones. Here, we are only 
concerned with the Fe-group composition. For the intermediate mass elements the main 
improvement compared to earlier calculations is the higher production of single nuclei like 
^^Sc and ^^V. Neutrino induced reactions become important for the nucleosynthesis of 
nuclei in the mass range A > 70. These nuclei are mainly produced in the zones close to 
the mass cut where the electron fraction strongly depends on the neutrino captures. 

4. Conclusions 

We have ensured explosions (of otherwise non-explosive models) by changing the neu- 
trino luminosity in two different ways. Detailed nucleosynthesis calculations based on a 
consistent treatment of all weak interactions show an electron fraction Ye ~ 0.5 for the 
innermost ejected zones as it is consistent with observations. Our calculations confirm the 
importance of consistent treatment of neutrino induced interactions for the nucleosynthe- 
sis of the innermost zones. In the Fe-group improvement was achieved for the isotopes 
58,6i,62]\jj^ 51 gj^^ ^^Sc. However, the nucleosynthesis calculations based on the innermost 

0. 03.0 remain to be extended to the outer layers to obtain a complete prediciton. 
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